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Abstract 

Influence of the plate surfaces roughness in precise ellipsometry experiments is 
studied. The realistic case of a Gaussian laser beam crossing a uniaxial platelet 
is considered. Expression for the transmittance is determined using the first order 
perturbation theory. In this frame, it is shown that interference takes place between 
the specular transmitted beam and the scattered field. This effect is due to the 
angular distribution of the Gaussian beam and is of first order in the roughness 
over wavelength ratio. As an application, a numerical simulation of the effects of 
quartz roughness surfaces at normal incidence is provided. The interference term is 
found to be strongly connected to the random nature of the surface roughness. 

1 Introduction 

The high-accuracy universal polarimeter (HAUP) jlj has proved to be a very useful in- 
strument to measure the crystal optical properties (see for instance 01111 El and references 
therein). The principle is simple and was introduced a long time ago (see [5] for an his- 
torical introduction): the light intensity measured after a rotating high quality polariser, 
a crystal plate (the sample) and a high quality rotating analyser, is fitted to a theoretical 
formula with several coefficients as free parameters where the delay due to birefringence 
and optical activity can be determined. 

The accuracy of this instrument has now reached the few 10~ 5 level and systematic 
errors contributing at this order of magnitude have been investigated jEHHIEJ- The con- 
clusion is that roughness is most likely one of the main source of systematic uncertainties. 

However, despite an extensive literature on surface roughness [HI QU|, no theoretical 
expression for the transmission of a Gaussian beam by an anisotropic rough platelet is 
available. It is the purpose of this article to provide this expression. We consistently 
take into account the Gaussian nature of the laser beam, the multiple reflection inside the 
plate and the roughness of both faces of the plate. To simplify the calculations we further 
restrict ourselves to uniaxial homogeneous crystals. As a result, we find that unlike plane 
waves, specular Gaussian beams are affected by the surfaces roughness, even in the first 
order perturbation theory. 
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The physical origin of this phenomenon is the angular distribution, or plane wave 
expansion, of Gaussian beams [TT] . Plane waves constituting a Gaussian beam having 
different wave vectors, a given plane wave can then be scattered in the specular direction 
of the other ones. The resulting interference pattern leads to an a priori non vanishing 
contribution of the scattered field in the specular region. To some extent, this phenomenon 
is thus related to the near-specular scattering by rough surfaces introduced in [T^j . 

Another aspect of realistic platelet surfaces is the interface parallelism default. De- 
pending on the wedge angle, this default can compete with roughness in the modifications 
of the transmitted beam polarisation. The nature of these effects is however different. 
Given the relative orientation of the two plate interfaces, the wedge effect is univocal 
whereas roughness, as it will be shown in this paper, is of random nature. It is then most 
likely that these two effects cannot compensate each other. In principle, the perturba- 
tive calculations reported in the present article holds for both effects. Nevertheless, the 
boundary matching method, applied a long time ago to isotropic wedges [13] , can be used 
to describe the wedge effect. We shall report this calulation in a future publication and 
restrict ourselves here on platelet roughness. 

This article is organised as follows. In section |21 we derive the theoretical expressions 
and numerical simulations of quartz plates are presented in section 



2 Formalism 

The choice of the theoretical formalism is driven by the properties of the crystal plate 
surfaces under study. Fortunately, an exhaustive experimental study on crystal surfaces 
has recently been published (Hj. Most of the high quality polished crystal surfaces used 
in optics have a profile surface correlation length of the order of the optical wavelength 
and a root mean square roughness of the order of a few angstrcem. It means that one can 
safely use a first order perturbation theory [T3] neglecting the local field effects [IE]. The 
more suitable formalism for our problem is the one introduced in and generalised to 
anisotropic overlayers in [TBI- However, in the latter reference, the anisotropy is treated 
perturbatively and only the reflection of plane waves is considered. We shall then extend 
this formalism to platelet 's transmission taking fully into account the plate anisotropy 
and treating perturbatively the plate roughness. 

In the following we tried to be concise, referring to f?l QB] for further details. The 
wave equation corresponding to the system represented in figure d is: 

V x V x E(r) = cjVoD(r) (1) 

with D(r) = £(z)E(r) and 

£{z) = e( -z + h (x,y) )e l + 



B[ z -h (x,y) ) - © ( z - a - h a (x, y) 

+e(z-a-h a (x,y))e l (2) 
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where 1 is the 3x3 identity matrix and G is the Heaviside function. For uniaxial media, 
it is useful to write [TH] 

e = ej_l + (ey - ejjcc 

with c T = {c x , c y , c z } the unit vector along the optical axis, cc a Dyad and e± = t? €q, 
eii = n^eo the ordinary and extraordinary components of the dielectric tensor. In equation 
(0), the two functions h (x,y) and h a {x,y) are the profiles of the two surfaces located at 
z = and z = a respectively. As usual [20], we assume that the two planes z = and 
z = a are defined such that the mean profiles vanish, i.e. < h >=< h a >= 0. 




Figure 1: Schematic view of the plane wave propagation in the anisotropic slab. For the 
sake of clarity, some of the inner reflected rays are represented by small arrows. The 
plane of incidence coincides with the plane yz. Symbols a± and b± correspond to the 
four possible propagation directions inside the medium. The vector basis {x, y, z} used 
throughout this article is also shown. The grey areas symbolise the scattered light due to 
surfaces roughness. 

The solution of equation (0) can be written E(r) = Eo(r) + E'(r) with Eo(r) given by 
the zero order wave equation 

V x V x E (r) = cuVoDo(r) (3) 
where D (r) = £ (z)Hj (t) and 

So(z) = e{-z)e l+ (e(z)-e(z-a])e + e(z-a)e l. (4) 
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To first order in uih/c [TH], one has £{z) = £q(z) + 5£(z) with 



8£{z) « h a (x,y)5(z - a) - h (x,y)S(z) 



(e± - e )I + (en - e ± )cc 



(5) 



with the Dirac distribution. 

To derive the differential equation for the first order scattered field E'(r), we introduce 
the Fourier transform 

E(K; z) = .F[E(r)] = — I E(r) exp(iK • R)d 2 R 

27T J 

where R = (r • x)x + (r • y)y and K = (k • x)x + (k • y)y with k the wave vector. Here, 
since we are considering Gaussian beams, no spatial length is introduced in the Fourier 
transformation. 

Taking the Fourier transform of equations JI} and (jHJ) and then subtracting them, one 
obtains fTB] 



^-zK + zJ^ ^-z'K • E'(K; z) 



dE' z (K;z) 

dz 



- -K 2 + 



92 
dz' 



rjE'QEC;*) =cu 2 fx B\K;z) 



(6) 

with D'(K; z) = D(K; z) — D (K; z). For perturbative stability, the wave equation must 
be written as a function of the continuous electric field components [T7], that is E' x , E' y 
and D' z . We shall do it separately for the s and p scattered waves as in ^7j. However, 
before providing the solutions we introduce |TH] the following useful vector function: 



F(r) = D(r) - £ (»E(r) = 5£(z)E(r) & F(r) = D'(r) - S (z)E'(r), 



(7) 



which gathers the infinitesimal contributions to the perturbated wave equation. To first 
order, one gets: 



FJt) 



FJt) 



1 



£zJz) 



£zz(z) 



£zz(z) 



5£ x Jz)DJ r) + 



£ zz (z) 8£ xx {z) - £ x Jz)b£ x Xz) 



EJ r) + 



£ zz (z)8£ xy (z) - £ yz (z)5£ xz (z) 



E y ( r] 



(8) 



5£ yz (z)D z ( r) + 



£ zz (z)5£ xy (z) - £ xz (z)5£ yz {z) 



EJ r) + 



£ zz (z)5£ yy (z) - £ yz (z)5£ yz (z) 



E y (r] 



(9) 



8£ z Jz)DJ r) + 



£ zz (z)5£ xz (z) - £ xz (z)5£ zz (z) 



EJ r) + 



£ zz (z)8£ yz (z) - £ yz (z)5£ zz (z) 



E y (r] 



(10) 



where £ij(z) and and 5£ij(z) are the components of the symmetric dielectric tensors 
of equations (0) and ©• In the leading order perturbation theory, one further set[T7] 
E x ( r) « E 0x ( r), E y ( r) ps E 0y ( r) and D z { r) « D 0z ( r) in equations JHE3U. 
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2.1 p scattered wave 



Projecting equation (0) onto z and utilising V • D = and equation 0, we obtain the p 
wave equation: 



d 1 dD' z (K;z) 
dzeo(z) dz 
A(z) 



cu 2 fi D' z (K;z)+ 
c z D' z (K;z) + e (z) 



K 2 



d 1 



_e (z)£ zz (z) 
D' Z (K; z) - S xz (z)E' x (K; z) - £ yz (z) t'jK: : 



c x E' x (K;z) + c v E , JK;z 



dz e (z) 



-iK-F(K;s) + i 



Hz) 

£zz(z) 



[K-c)c z F z (K;z) 



+ 



where we introduced e (z) 



e(z)-e{z-a) 



e(-z)+e(z-a) 



e(z)-e{z-a) 



e± and A (z) 



e || _ e _i_) such that equation (0J reads £0(2) = e o(^)I + A(z)cc. 

Solutions of equation (fTTj) are obtained using the Green's functions[T71 [TH]. There 
exist, a priori, nine Green's functions and thanks to the Dirac distributions appearing in 
equation ©, they must only be determined for z' = and for z' > a (we remind that we 
are interested by the solution in the region z ^> a). Furthermore, for z 3> a and z < 0, 
all terms of equation (fTTj) is front of the field components E' X (K; z) and E' y (K; z) vanish. 
Hence, the wave equation being expressed as function of the continuous field components, 
only one non zero Green's function G p (K; z, z') exists [T71 ITS] in the two relevant regions 
z > a', z' < and z > z', z' > a. The solution of equation (JTJJ) therefore reads: 



D'JK:z) 



G p (K;z,z') 



d 



dz' e (z>) 



—iK ■ F(K: z') + i 



S zz (z') 



'K-c)c z FJK:z') 



for z > a, where the Green's function is given by [21] 



G p {K;z,z') 



± ^<(K; z)E>(K; z')Q(z' - z) + £>(K; z)£<(K; - a 7 )) 



(12) 



(13) 



with 



dE>(K;z) 
W = E<(K:z) p( 1 



E>(K;z 



dE<(K;z) 



dz p v ' ' dz 

according to a theorem that can be found in [22] ■ Here Ep(K;z) and Ep(K]z) are 
the two independent plane- wave solutions of the unperturbated equation (J5J): E p > (K.;z) 
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corresponds to a wave coming from z — > — oo and Ep(K;z) to a wave coming from 
z — > +00. They are thus defined by the following boundary conditions: 

lim i?r"(K; 2;) oc exp(— ik z z) 



z— >+oo 

lim £^(K; z) oc exp(i/c z z) . 



with k z = + (k 2 -K 2 ) 1 ' 2 . 

Integrating by part the first term in the integral of equation (|T2j) and using equations 
PEED}, one gets 



£J(K;«) =W -zKcA n 1 + 



1 



c 2 ,D? (K) 



— 1- 



n 2 {n 2 + A n c 2 z ) 2 
., ( K x££ (K; 2) + K • y££ (K; z 



ni - 1 / - 



77. 



A n K • c 



x ( M'(K)c x + H%\K)cy 
1 



+ K 



1 + 



< - 1 + A n c z z 



D z {K;z) 



+ c z A r 



c x E p x (K; z) + CyE%{K; z) 



(14) 



for z > a and with A n = n 2 — n 2 and K = (K x x + K y y) / K. To obtain this expression, 
we used the definition Ep(K; z) = k / 'K~E' z (Jt; z) with K = y when K = ^H] (see figure 
^for the definition of the reference axes). To shorten equation (JUjl we also introduced 



D Z (K; z) =G p (K;z } a)F 

G p (K;z,0)F 
dG p (K; z, z') 



D 0z (R\ a) 



hJR) 



A>*(R; 0) 



A 

MR) 



Dl (K:z) 



dz' 

dG p (K;z, z') 



dz' 



T 

1 



Aiz(R; a 



z'=0 



D 0z (R; 0) 



h a (R) 
A 

ho(R) 



(15) 



(16) 



where A is the laser wavelength. Identical expressions hold for E X (K), Ey(K), E x (K) 



and Ey (K) 



To derive equation ()14j) we assumed that 



f(z)5(z)dz = - 



lim f(z) + lim f(z) 

2^0 + z~>0~ 



(17) 



where f(z) is a discontinuous function, but with a finite jump. Although this expression 
is not mathematically justified as stated in [22], it can however be used by considering 
that the Heaviside functions of equations (I2I4|1 are given by the limit 

G(z) = Km[l + tanh(>/C)]/2 . 

This choice is justified by the freedom existing in the determination of the dielectric 
tensor at z = [21]. It is to mention that equation (fTTj) leads to a disagreement with 
the boundary matching method for isotropic-isotropic interfaces in the case of oblique 
incidence. Another prescription was proposed in [23] to avoid this discrepancy. But, as 
mentioned in J2E], no general proof was provided in (2H|. There is then no reason for this 
particular prescription to work also for isotropic-anisotropic interfaces. In addition, since 
we are going to restrict ourselves to normal incidence, we choose to use the more intuitive 
and symmetric prescription of equation (JT7j) for our calculations. 



2.2 s scattered wave 

Following the lines of the previous section, we get the s wave equation: 



rj2 

K 2 -u 2 t , e (z)- — )E' s (K;z 



dz< 



A(*) 



-s ■ c e {z 



c x E' x (K;z) + c y E' y (K;z 



+ c z D' z (K;z 



A(z) 



a//i ( F(K; z) ■ s - ■ sc 2 F 2 (K; z) 



;is) 



with E' S (K; z) — s ■ E'(K; z) and § = {—Kyk. + K x y)/K. The solution is given by 

^(K; z) = u 2 v J°° G S (K; z, z') ^F(K; z') ■ s - ^^y£ ■ sc 2 F 2 (K; dz' , 

with G S (K; z, z') the s wave Green's function for which an expression similar to equation 
(I13|) holds. After integration, one finds 



E'(K-z) 



2tt 5 



A 



y - 1; 

| An 

n 2 Q + A n c 2 z 
+ A n s -ell 



s x E x (K; z) + s y Ey(K; z 
A 2 c 2 

n z 



s • c \ n 



n 2 Q + A n c 2 z 

1 



c x E x (K; z) + c y E y (K; z) 



(n 2 + A n c 2 z ) 2 



c z D z (K;z) 



(19) 



where E%.(K;z), Ey(K;z) and D S Z (K; z) are obtained by substituting G P (K; z, z') by 
G S (K; z, z') in equations (jl5H16|) . 
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2.3 Transmitted intensity 



Anticipating the numerical studies of section 01 we shall now consider a Gaussian beam 
at normal incidence coming from the region z < 0. Expressions for the electric field at 
z = and z = a read: 



E (K;0) = E i (K;0) + E r (K;0) 
E (K;a) = E t (K; a) 



with|27j 



E,(K;0) = -^=exp(-^- )M 3x3 e (20) 
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E t (K;a) = -^exp( -^f- )QM>Q T M 3x3 e (22) 



E r (K;0) = -^exp( ^- ]OM>fl J M 3x3 e„ (21) 

V 27T 



u> ( wlK 2 
/2tt 

where we chose the beam waist position at z = and where |2*7]: eo is the electric vector 
describing the polarisation of the Gaussian beam centre (i.e. K = 0), M 3x3 is a 3 x 3 
matrix describing the polarisation of the plane waves constituting the Gaussian beam 28 , 
M> and Mf are the Jones matrices describing the reflection and transmission by the 
uniaxial parallel plate (the upper script > indicates that these matrices correspond to an 
incident wave coming from z < 0). The Jones matrices take into account the multiple 
reflections inside the platelet. They are determined [2Hj in the basis {§, p, k} and then 
transformed to the basis {x, y, z} thanks to the transfer matrix Q. 
The Green's functions are given by 



G m (K; z, z' 



(Zih)- 1 ^exp(ik z [z -z'}) + M< mm exp(ik z [z + z']) ) ; z' > a 
X2ik z )- 1 M< m exp(ik g [z-z']);z' < 



with m — 1,2 for s and p waves respectively. It is to mention that M< = M> and 
Mf- = M> when the optical axis is in the plane of interface, i.e. when c z = 0. 

From the above expressions, one can compute equations (|15H16|) and then the s and p 
scattered fields. In doing so, the following kind of Fourier transform is to be evaluated: 



E 0x (R; a)h a (R) 



j J E 0x (K';a)h a (K-K f )d 2 K', (23) 



where, because < h a >= 0, the Fourier transform of the surface profile h a (K — K') 
vanishes when K = K'. However, and this is one of the major points of this article, 
since £oa;(K;a) oc exp(— WqK 2 /A) at normal incidence, then ^ 7 [£ l o a; (R; a)/i (R)] does not 
necessarily vanish when K = as it is the case for a single plane wave. Consequently, 
the specular transmitted beam receives a non vanishing contribution from the scattered 
field, even in the first order perturbation theory. If the Gaussian beam is viewed as 
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a superposition of plane waves [TT], then this phenomenon is then due to the angular 
distribution of the plane waves. However, this contribution depending on h a (K — K 7 ), a 
realistic simulation of the surface roughness is needed to evaluate the integral of equation 
(I23J) . This is the subject of the next section. 

To exhibit this contribution, let us assume that a Wollaston prism is located after the 
anisotropic plate and that its axes correspond to the x and y directions. Writing the 
scattered electric field in the basis {x, y, z}, 

E' X (K; z) = fiu££(K; z) + n 12 E' p (K; z) 
E' y (K; z) = n 2l E' s (K; z) + n 22 E' p (K; z) 

one obtains 



d 2 K 





x-E t (K;z)| 2 d 2 K + J J |^(K;^)| 2 d 2 K+ 
x ■ E t (K; z) (e' x {K; z)) + U ■ E t (K; z)\ E' X (K; z) 




|y-E t (K;^)| 2 d 2 K + 




|£'(K^)| 2 d 2 K+ 



J J y-E t (K;z)^(K;^ + ^y • E t (K; z)^j E' y (K; z) 



d 2 K 



for the two intensities I x and I y measured after the Wollaston. Writing I x = ll +/| +/| , 
and I y in the same way, with 



m 

X 

J[l] 

X 

j[2] 




|x-E t (K;^)| 2 d 2 K 





x-E t (K;z)f^(K;z; 
E'(K-z)\ 2 d 2 K 



x.E t (K;z) E' x (K;z 



d 2 K 



(24) 
(25) 
(26) 



one sees that the specular-scattered interference term Ix is of first order in a/ < h 2 >/X. 
As wq — > oo, the Fourier transforms of the electric field components lead to Dirac distri- 
butions and the usual plane wave result is recovered. This interference term is therefore 
expected to be dependent on the laser beam waist. 

Let us finally note that equation (J2lj|) does not completely describe the second order 
contribution I x in the specular region, the interference between the specular and the 
second order scattered field being omitted in our calculation. 



3 Numerical simulations 

To estimate the specular-scattered interference term, a laser beam crossing a quartz 
platelet at normal incidence is now considered. The incident electric vector is fixed to 
e = x, i.e. perpendicular to the plane of incidence of the Gaussian beam's centre, and 
the intensities recorded after a Wollaston prism are calculated as in the previous section. 
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3.1 Numerical ingredients and input parameters 

As described in the previous section, the specular-scattered interference term can only 
be evaluated if a simulation of the surface roughness is provided. The profiles of the 
quartz surfaces ho and h a are thus randomly generated independently. Then the scattered 
fields are computed for various orientations of the optical axis and the two intensities 
of equations (j24H26J) are determined. The procedure is repeated in order to obtain a 
statistical distribution of the intensities. The numerical integrations of equations (j24H26|) 
are performed in the domain arctan(K/k) < 1°, which roughly corresponds to the angular 
acceptance of a 1 inch diameter optical lens located at 1 m from the plate. In practice 
this limit does not affect the value of the interference term but only the second order 
contribution. 

Random profile generations are performed as follows. First, the height distributions 
are determined according to a Gaussian distribution of mean value zero and root mean 
square a = 6 A. This is a typical value for a high quality polished quartz plate |3D~] . The 
heights are stored in a grid {x, y} of length L x L containing 2 n x 2 n nodes. The value 
of n is limited by the computer memory capacity, n — 11 in our case. Next, the Fourier 
transform is computed and the profile heights are weighted by the square root of the 
two dimensional spectral density function [20] (PSD2). Though we are going to consider 
an isotropic roughness distribution, we shall not use the reduced one dimensional radial 
PSD1 20JEI]. In doing so we fully account for the random nature of the surface roughness 
in the evaluation of the specular-scattered interference term (the use of a PSD1 would 
leads to a smaller dispersion of our numerical results). 

The result of jSHj for the PSD2 is used: 

PSD2(K) = { V ; ' - mm 27 

[0;K<K min 

where the cut-off spatial frequency K min has been introduced to account for the surface 
profile property < h >= => PSD2(0) = [22]. The correlation length I is of the order 
of[30 j 1/im and K m i n is smaller than 10 _3 /im _1 [14] . The parameter 1/K m i n acts as a 
spatial frequency threshold for the laser radius: roughly speaking, for wo > \/2/ K min the 
Gaussian beam behaves as a plane-wave and the specular-scattered interference term van- 
ishes. Nevertheless, the present values of K min are limited by the acceptance of the surface 
profile measurements: values as small as 10 _4 /im^ 1 for Si substrates [32J and 10 -5 /zm _1 
for K min for Si wafers [33J have been reported (notice that these numbers lead to a laser 
radius threshold greater than 15 mm). Finally, the inverse Fourier transform is computed 
leading to a 'coloured' random surface. 

The grid parameters are determined by the correlation length I and the laser waist 
wo. The distance between two nods of the grid {x, y} must be at least twice smaller than 
I and the distance between two nods of the grid {K x , K y } ten times smaller than \/2/wq 
approximately. This leads to the following 'experimental' parameters: Wq = 100/im, 
I = 1.6/im and L = 8wq. As for the laser wavelength we choose A = 0.6328 /im and 
n = 1.542637 and n e = 1.551646. The ratio a/A is therefore of the order of 10~ 3 in our 
numerical examples. 
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Using the numbers given above, we have written a computer program to estimate the 
specular-scattered interference term. Calculations of the unperturbated electric fields and 
of the Green's functions are performed according to A fast Fourier transform (FFT) 
algorithm is used for the simulation of the surface roughness and for the integrals of the 
kind of equation (|2H|) . The numerical precision for the unperturbated intensities is of the 
order of 10 _6 j2Z]. The specular-scattered interference term is thus known to an accuracy 
better than a/A x 1(T 6 « 1(T 9 (with our grid size, the FFTs do not reduce this accuracy). 
However, the CPU-time is quite sizable: with a SPECfp2000 1288 computer, the random 
generation of the surfaces takes 22 s and next the computation of the intensities for one 
given orientation of the optical axis takes 183 s. 

3.2 Numerical results 

We first consider a quartz plate thickness a = 0.720 mm with the optical axis located in 
the plane of interface (9 C = arccos(c 2 ) = 7r/2), that is a tenth order quarter- wave plate. 
1^ and I y °^ are shown is figureElas function of the optical axis azimuth C = arctan(c 3/ /c a; ). 




0.5 1 1.5 2 2.5 3 

<\> c I rad 



Figure 2: Intensity of the specular beams as function of the optical axis azimuth for 
various quartz plates: a tenth order quarter wave plate (full dots), 0.562 mm thickness 
(open dots), 5 mm thickness (diamonds) and a tenth order plate thickness but with the 
optical axis polar angle fixed to 9 C = 7r/4 (stars). The upper set of curves represents the 
intensities ll and the lower set of curves I y °K 
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The results for the first order contributions I x and I y are shown in figures X\\\\ 




§ c I rad 

Figure 3: Interference between the specular and scattered transmitted fields ll as func- 
tion of the optical axis azimuth. The plate is a tenth order quarter wave plate and the 
beam waist is wq = 100 /iin. Different symbols correspond to different random generated 
surface profiles. 

Each curve of these plots corresponds to different surface profiles. Considering one 
given profile, one can notice that: the size of the specular-scattered interference strongly 
depends on the surface profile and can reach the per mill level of the zero order contribu- 
tion, its sign changes with C and its shape is not regular with C . The change of sign is 
expected since the intensity averaged over a large number of profiles obviously vanishes. 
The erratic shape is also expected since the fields change with C and so do the Fourier 
transforms as the one of equation (J2*3*|) . 

The second order contribution (calculated from equation (j2EJ)) is six order of magnitude 
smaller that the zero order contribution. However, we do not show any results since our 
second order calculation is not complete concerning the specular angular range. 

Large differences are indeed observed when the plate thickness is changed. The 
specular-scattered interference contributions are computed for a = 0.562 mm as in jH] (i.e. 
(8 + 10 -3 ) x 27r retardation plate with our choice for the optical indices) and a = 5 mm 
as in |7j (i.e. (71 + 0.18) x 2n retardation plate), and still with 6 C = tt/2. They are 
compared to the values obtained with the tenth order quarter wave plate and the same 
surface profiles. 
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Figure 4: Interference between the specular and scattered transmitted fields I y as func- 
tion of the optical axis azimuth. The plate is a tenth order quarter wave plate and the 
beam waist is Wq = 100 /iin. Different symbols correspond to different random generated 
surface profiles. 

The results are presented in figures loIHl 1% and Iy scale with j| ' and -4°' (see figure |2J). 
In particular, the oscillations of l| are dumped when j| ' gets flat (i.e. for the almost 
zero retardation plate a = 0.562 mm). 

To investigate the dependence of equations ()19|14j) on the optical axis polar angle 9 C , 
the calculations were performed fixing 8 C = ir/4 for the plate thickness a = 0.720 mm. 
Here again the variations are noticeable (see figures 03 and [HJ) . 

Looking at figures H] and one can remark that I y tends to be of opposite sign in 
the regions < C < 7r/2 and 7r/2 < <p c < n. But this is not a general rule as it seems 
to come out from experiments [HI Ej . One can also observe two fix points at <p c = and 
7r/2 on figures |U and El Iy being the interference between the scattered field and the 
zero order field, these fix points correspond to the zeros of the zero order field (see if^ 
on figure |2J). This is not the case for the second order term of equation (|26[) which is 
of the order of 10 -6 and therefore dominates around C = 0,7r/2 (here the missing term 
of equation (j2*o^l is not relevant since it describes the interference between the specular 
and the second order scattered fields). However, the dispersion of around zero for 
0c = 0,7r/2 (not visible on these figures) defers very slightly from zero, it is of the order 
of 10~ 10 for wq = 100 fim and ~ 10~ 8 for wq = 25 /im. This is a cross-polarisation effect, 
i.e. this is due to the matrix M^ x ^ in equations (j20H22j) . 
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Figure 5: Interference between the specular and scattered transmitted fields I x as func- 
tion of the optical axis azimuth for: a tenth order quarter wave plate (full dots), 0.562 mm 
thickness (open dots), 5 mm thickness plate (diamonds) and a tenth order plate thickness 
but with the optical axis polar angle fixed to 6 C = tt/A (stars). The surface profiles are 
the same for the four curves. 

The numerical results presented here are rather independent of the choice for the 
PSD2 provided a quartz plate of optical grade is considered. It is indeed experimentally 
demonstrated [TJJ EDI Ell EH ESI that the PSDs of optical element's surface have an 
inverse-power-low (or Fractal-like) behaviour. Therefore, various smooth mathematical 
representations of the PSD (see [HH EH] for examples) are reducible to equation (J27)l as it 
is justified in p?T] . 

As a concluding remark, we point out that three important dimensionless parameters 
LUo/l, a/X and l/X have been encountered in our calculations (ujq/X describes the cross- 
polarisation effects discussed above and is therefore not relevant here). 

As mentioned in section El the validity of the perturbation treatment depends on a/X 
and l/X. Using the first order perturbation theory, the severe conditions <x/A <C 1 and 
l/X ~ 1 must hold[T3]. They are fortunately fulfilled by optical grade quartz plates. To 
study the influence of the correlation length I, we changed the value of / to 0.7/iin and 2//m 
and we observed no significant qualitative differences with respect to the results described 
above. 



14 



-0.6 

0.5 1 1.5 2 2.5 3 

<|> / rad 



Figure 6: Same as figure but for I y . 

As for the last dimensionless parameter ooq/1, we already mentioned that when ujq/1 — > 
oo the usual result for plane waves is recovered (i.e. the specular-scattered interference 
term vanishes) although, with regard to the cut-off parameter K m i n , this limit seems to be 
idealistic for a finite size quartz plate. To get an idea of the influence of wq, we increased it 
to 200/mi and here again, no significant differences were observed. Much larger values for 
wo could not be tried, keeping a reasonable correlation length, because of the computer 
memory limitation. Finally let us mention that the other limit Uq <C / corresponds to the 
scattering by gratings jHH] • I n this limit the specular-scattered interference term vanishes 
since the diffusion occurs at large angle with respect to the specular beam direction. 

4 Conclusion 

We have computed, in the leading order perturbation theory, the effect of surfaces rough- 
ness on uniaxial platelets transmittance. Taking into account the Gaussian nature of 
laser beams we showed that the interference between the specular and scattered fields 
contributes to the intensity measurement performed in the specular region. 

This contribution is of first order in the ratio of the root mean square roughness over 
the laser wavelength cr/X. It depends strongly on the plate surfaces profiles and on the 
crystal optical properties, orientation of the optical axis, thickness and optical indices 
(i.e. temperature). It is therefore useless to implement the roughness calculation in a 



15 



HAUP type of fitting procedure (in addition, the numerical calculation are computer 
time consuming). 

In view of our numerical results, it is most likely that simple overlayer models [E] can- 
not describe accurately the dynamical properties of our main formula equations (|14|19j) . 
Nevertheless, we point out that, because of the random property of the specular-scattered 
interference term, a simple way to avoid it is to perform a series of measurements at 
different locations on the plate and then to average the results. Although this proce- 
dure would increase the uncertainty on the determination of crystal optical parameters, 
it should however decrease the systematic bias. The determination of the plate thickness 
in situ could be done by varying the laser incident angle (i.e. by tilting the plate) [37] . 
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